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Abstract. We investigate the possibility and stability of bandferromagnetism in the 
single-band Hubbard model. This model poses a highly non-trivial many-body prob- 
lem the general solution of which has not been found up to now. Approximations 
Q_) ■ are still unavoidable. Starting from a simple two-pole ansatz for the spectral density 

our approach is systematically improved by focusing on the influence of quasiparticle 
damping and the correct weak-and strong coupling behaviour. The compatibility of 
the different aproximative steps with decisive moment sum rules is analysed and the 
— < ' importance of a spin-dependent band shift mediated by higher correlation functions 

is worked out. Results are presented in terms of temperature- and band occupation- 
dependent quasiparticle densities of states and band structures as well as spontaneous 
>-rt \ magnetisations, susceptibilities and Curie temperatures for varying electron densities 

<*"] ■ and coupling strengths. Comparison is made to numerically essentially exact Quantum 

O ' Monte Carlo calculations recently done by other authors using dynamical mean field 

O , theory for infinite-dimensional lattices. The main conclusion will be that the Hubbard 

model provides a qualitatively correct description of bandferromagnetism if quasipar- 
ticle damping and selfconsistent spin-dependent bandshifts are properly taken into 
account. 

m 

CN . 1 Introduction 

t> 

Ferromagnetism is bound to the existence of permanent magnetic moments. If 
these belong to itinerant electrons in a partially filled conduction band one speaks 
of bandferromagnetism JjJ . Archetypical representatives are the classical 3d fcr- 
romagnets Fe, Co, Ni. The microscopic interpretation of bandferromagnetism 
is one of the most interesting and most complicated many-particle problems in 
solid state physics. The simple but fairly successful "Stoner criterion" 

g ! Uqo(E f ) > 1 (1) 

(U: intraatomic Coulomb interaction, qq: free Bloch-density of states (B-DOS), 
Ep- Fermi energy) defines as a minimum set of ingredients for a theoretical 
model to describe bandferromagnetism the Pauli-principle, the kinetic energy, 
jJJ ■ the Coulomb interaction (strong and strongly screened), and the lattice struc- 

ture. This minimum set is realized in the Hubbard-Hamiltonian (2) for correlated 
electrons on a lattice in a non-degenerate energy band, 
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H = J2( Ti 3 - f iS v) C t C 3° +2 U ^2 n i° n i-° ( 2 ) 
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(/z: chemical potential). The Coulomb interaction is restricted to its intraatomic 
part only (rii a — c^Cia), while the kinetic energy contains the hopping inte- 
grals Tij being strongly related to the lattice structure. The Pauli principle is 
guaranteed by the formalism of second quantization. cf a {ci a ) is the creation 
(annihilation) operator for an electron with spin f or 1 at lattice site Ri. 

The physics of the Hubbard-model is decisively determined by the Coulomb 
coupling U/W (W: Bloch-bandwidth), the lattice structure and the band occu- 
pation n = jj J2a < n a > (0 < n < 2; N: number of lattice sites). 

In spite of its simple structure the Hamiltonian (2) provokes a rather sophisti- 
cated many-body problem, that could be solved in the past only for a few special 
cases. In general approximations are unavoidable. So it was not fully clarified 
until recently whether or not fcrromagnetism is possible in the Hubbard-model 
for finite U, finite temperatures T and band occupations away from half filling 
n=l. Nagaoka showed p[ for the special case U — oo and electron numbers 
N e = N ± 1 on an sc or bec lattice and N e = N + 1 on an fee lattice, respec- 
tively, that the fully spin-polarized particle system represents the ground state. 
However, the (numerical) proof of finite-temperature ferromagnetism in an ex- 
tended parameter region was only recently given by Ulmke H for an infinite 
dimensional (d = oo) fec-type lattice by use of Quantum Monte Carlo (QMC) 
calculations in connection with Dynamical Mean Field Theory (DMFT) JUfJ . On 
the other hand, the certainty that ferromagnetism exists in the Hubbard model 
does not at all mean that the phenomenon itself is understood. What is the phys- 
ical mechanism enforcing ferromagnetic spin order of itinerant electrons? This 
question can be answered, if at all, better by analytical approaches to the highly 
complicated electron correlation problem than by purely numerical evaluations. 
Starting from some basic features we shall construct in the following chapters 
three different theories for the Hubbard model to contribute to an answer of 
the above question. The theories are constructed in such a way that each differs 
from the preceding one by eliminating an obvious shortcoming. We start from a 
simple spectral density approach (SDA) which suggests itself by some rigorous 
strong coupling features. However, it suffers from a complete neglect of quasi- 
particle damping and a wrong weak-coupling behaviour. By use of a modified 
alloy analogy (MAA) the advantages of the SDA are retained but quasiparticlc 
damping is included. Comparison of the SDA- and the MAA-results helps us 
to recognize the influence of quasiparticle damping on magnetic stability. The 
weak coupling behaviour, however, remains insufficient. From a modified per- 
turbation theory (MPT) which correctly reproduces the weak- as well as the 
strong-coupling limit, we learn by comparison with SDA and MAA how decisive 
the correct weak coupling (low energy) description is for a theoretical approach 
that aims at the strong-coupling phenomenon ferromagnetism. 



2 Hubbard model basic features 

It is commonly accepted that ferromagnetism is a strong-coupling phenomenon. 
A theoretical approach should therefore be reliable first of all in the U ^> W 
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regime. Practically all information we are interested in can be derived from the 
retarded single-electron Green's function H 

C k ,-Ai , , -i / #exp(iEt)<[<*«r(*), <&(«)]+) ( 3 ) 
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Cka 



—= Y^ c ia exp(-ifcHi) (4) 



[.., ••]+(-) means the anticommutator (commutator) and (..) the grandcanonical 
average. A; is a wavevector from the first Brillouin zone. Of equivalent importance 
is the single-electron spectral density being directly related to the bare line shape 
of an angle- and spin- resolved (direct or inverse) photoemission experiment: 

S k a{E) = --ImG ka (E + i0+) (5) 

According to the pioneering work of Harris and Lange M we know that in the 
strong coupling limit Ska (E) is built up by two main peaks near To and To + U 
where To is the centre of gravity of the free Bloch band, and additional satellite 
peaks at T + pU {p = — 1, —2, . . . ;p = 2,3, . . .) (s. Fig.l). The spectral weights 
(ex peak area) of the satellites, however, decrease rapidly with increasing distance 
from the main peaks. Already the next neighbours (p = — 1, +2) acquire only a 
weight of order (^j- ) 4 being therefore negligible in the strong-coupling regime. 
So we can assume for U ^> W a two-peak structure of the spectral density. The 
exact shapes of the peaks are not known, but the centres of gravity are, 

T la (fc) = T + (1 - (n- a ))(e(k) - T ) + (n_ ff )B fc _ ff + 0( — ) 4 (6) 

W A 

T 2 a{k) =T Q + U+ (n_ CT )(e(fe) - T ) + (1 - (n_ a ))B k _ a + 0( — ) 4 (7) 
as well as their spectral weights: 

W 

a la (k) = (1 - (n_ CT » + 0( — ) - 1 - a 2a (k) (8) 

At least formally the centres of gravity carry a spin-dependence which may give 
rise to an additional exchange splitting of the two main peaks as the fundamental 
precondition for ferromagnetism. In the following it is demonstrated that such 
a spin asymmetry is mainly due to the "band correction" -Bfc_ CT : 

-Bfc-cr = B- a + Fk- a (9) 

The local term B^ a ("band-shift") can be interpreted as a correlated electron 
hopping: 



(n_ CT )(l - (n- ff ))B- CT = -^T ij (c l t CT c J _ ff (2n iCT - 1)) (10) 
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Fig. 1. Single electron spectral density of the Hubbard model in the strong coupling 
regime (U 3> W) as function of energy. 



For U 3> W and less than half-filled bands double occupancies are very unlikely 
so that the second term in the bracket dominates. That means that the shift 
of the (T-spectrum is correlated with the negative kinetic energy of the (— <r)- 
electrons. For further evaluations it will turn out to be decisive that the band- 
shift -B_ CT can exactly be expressed by the single-electron spectral-density [||^] 

(n_ CT )(l-(n_ (7 ))B_ (7 = 

1 f°° 2 

— J2(e(k)-T Q )J dEf.(E)(-(E-e(k))-l))S h . a {E-n) (11) 

f-(E) = (exp(3(E — /i) + 1) _1 is the Fermi function. B_ a obviously disappears 
in the zero-bandwidth limit(e(fc) — ► ToVfc). 

The k-dependent part of the band correction (9) is built up by density- 
density, double hopping and spinflip correlation terms 

{n_ a )(l-{n_ a ))F k _ a = 



N 



+ ( c j-<j c ja c i-° c i<?) + (Cjo-Cj-vCi-aCiv)) 



(12) 



It vanishes in the zero-bandwidth limit and has no direct influence on the centres 
of gravity T\^, a of the quasiparticle subbands (Hubbard bands) that emerge from 
the two spectral density peaks (Fig. 1): 



k 

T 2CT - -j= J2 T 2a (k) =T + U+(1- (n_ CT »B_ CT 



(13) 

(14) 
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Fk-a may, however, lead to a spin-dependent bandwidth correction competing 
then with the other band narrowing terms in (6) and (7), respectively. The im- 
portance of Ffc_ (T with respect to ferromagnetism shall be discussed in Sect. III. 
Unfortunately, it cannot be expressed by the spectral density as the band shift 
B- a by eq. (11). A determination of F k - a therefore requires further approxima- 
tions. 

If -Bfc-er has indeed such vital implications a theoretical approach should 
handle this term with special care. That can effectively be controlled by the 
spectral moments of the spectral density, 

oo 

M kJ = \ J dEE n S k „(E) n = 0, 1,2,3, . . . (15) 

— OO 

which can be calculated independently of Ska{E) via 

M tJ = (l(rt§- t TcUt),cL(t')} + )t=ti (16) 

In practice, however, the moments are useful only for low order n because with 
increasing n eq.(16) produces higher expectation values which are usually un- 
known and not expressible by the spectral density, either. The important band 
correction B k _,j(E) first appears in M^J ||pl. 

<J = 1 (17) 

M^ = e(k) + U(n_ a ) (18) 

M kJ = £ 2 ( fc ) + 2s(k)U(n- a ) + U 2 (n^) (19) 

M k J = i 3 (k) + 3e 2 (k)U(n_ a ) + e(k)U 2 (n_ a )(2 + (n_ a )) (20) 
+t/ 3 (n_ CT ) + U 2 (n- a )(l - (n- a ))(B k _ a + T - fi) 

(i(k) — e(k) — /i). A necessary condition for a theoretical approach to be con- 
sistent with the strong-coupling results (6), (7), (8) is that the first four moments 
MuL , n = 0, 1, 2, 3 are correctly reproduced. The condition becomes even suffi- 
cient when additionally the zero-bandwidth limit |l(| is fullfilled. One can ele- 
gantly check the strong coupling consistency by use of the high-energy expansion 
of the single-electron Green's function 

GUE)= (dE'^^n^^ (21) 
«=o 

which can also be used in the Dyson equation, 

GkAE) = G k 0) (E) + G k 0) (E)^S ka (E)G ka (E), (22) 

where 

Gf{E) = h(E + ^-e(k))- 1 (23) 
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is the U = Green's function, to yield a respective expansion for the selfenergy: 

S UE) = £ ^ ( 24 ) 

m=0 

The coefficients Cj^ — f m (M^J , . . . , M^™ ) are simple functions of the mo- 
ments up to order m + 1: 

CJ®=Mi»-(e{k)-ri (25) 

C%=M%-(Miy (26) 

C ( £ = M<% - 2M&>M£? + (M«) 3 (27) 



For a given theoretical approach one can easily control whether or not the lead- 
ing term in the expansions (21)and (24), respectively, match with the rigorously 
calculated moments. However, even the opposite procedure may be useful, infer- 
ing approximately the overall energy dependence of the Green's function and the 
selfenergy, respectively, from the first few exactly calculated terms. A proposal 
how this can be done is exemplified in the next Section. 

Rigorous statements are of course also possible in the weak coupling regime 
(U <C W). As mentioned in Sect. I it remains a challenging problem to find out 
how decisive weak coupling and low-energy properties are for a theoretical ap- 
proach to correctly describe the strong-coupling phenomenon "ferromagnetism" . 
In the "diagram-language" the Dyson equation (22) reads 



■ f*\ . 



Gfa(E) = G'°\E) + Gf{E)^E k „(E)G k „{E) 

Using standard perturbation theory [0,0,0 one has to sum up for the selfen- 
ergy Sk a (E) all dressed skeleton diagrams. A "skeleton diagram" is a selfenergy 
diagram that does not contain any selfenergy insertion in its propagators. If in 
addition the propagators are "full" Green's functions Gktr{E) then the skeleton 
is "dressed". That means for the Hubbard model up to second order: 
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Conventional diagram rules [p 1 1 , 12] yield as selfconsistent second order pertur- 
bation theory in-^: 

Z k *{E) = U(n- a ) + U 2 A ka (E) + 0(U 3 ) (28) 

The linear term is the Hartree-Fock (Stoner) part while the more complicated 
second order contribution is in general non-local, complex, and energy depen- 
dent: 

A k a(E) = 

Wr^j J J d£ld£2d£3 E + v-e 1+ e 2 -e 3+ m 

x(/_(ei)/_(-e a )/_(e 8 ) + /-(-ei)/-(ea)/-(-e3)) (29) 



S is the full spectral density which has to be determined self-consistently Several 
alternatives are thinkable which are correct and equivalent up to order U 2 . One 
of these is to replace the "full" spectral density in (29) by its Hartree-Fock 
version jl|Q 

S£(E)=S^(E-U{n-„)U) (30) 

To avoid an ambiguity we require that the Hartree-Fock particle densities are 
the same as those in the full calculation: 

(n CT )« = (n a ) (31) 

that can be regulated by considering the chemical potential fi^ in (30) as a 
proper fit parameter. There are of course other possibilities to avoid the men- 
tioned ambiguity p5| . 

Several very important conclusions can be drawn from (28), e.g. for the low- 
energy behaviour of the selfenergy Q at T — (a, /3, 7, real), 

S ka (E) = a ka +(3 ka E + i 7krT E 2 + ... (32) 

That means that the imaginary part of S disappears at the Fermi edge (E = 0) 
indicating quasiparticles with infinite lifetimes. At low but finite temperatures 
holds 

ImS ka {0) ex T 2 (T — ► 0) (33) 

(32) and (33) are necessary to guarantee the correct Fermi liquid behaviour. 



3 Analytical Approaches 

We can now use the basic features of the last section to construct and control 
three different analytical approaches. The main goal is to work out the essentials 
for ferromagnetism in the Hubbard model by a critical comparison of these 
theories. 
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3.1 Spin-dependent band shift 

Having in mind the two-peak structure of the spectral density in the strong- 
coupling regime (Fig.l) we can construct a simple "spectral density approach" 
(SDA). If we initially assume that quasiparticle damping, responsible for the 
finite peak widths, does not play a dominant role for magnetic properties then 
the following two-pole ansatz appears to be plausible: 

2 

S ka {E) = nJ2aja{k)6(E + n-E j(T (k)) (34) 

i=i 

The quasiparticle energies Ej a (k) and their spectral weigths aj a (k) are easily 
fixed by equating the first four moments. The resulting sclfenergy has a remark- 
able structure (|,g,|l^,^8|, in particular the important "band correction" Bk- a 
is involved: 

S ka V E ) = U{n- a )-— - -— (35) 

E + n — T — B k - a - (7(1 - {n-a}) 

According to (10) and (12) B k - a vanishes in the zero-band-width limit (W — ► 



0). £%» A (E) is then the exact W = selfenergy [pM7|. Additionally the SDA 



fullfills by construction the first four moments so that the strong-coupling be- 
haviour must be correct according to (6)-(8).The first three terms of the high- 
energy expansion of U^ (E) agree exactly with (25)-(27). 

We note in passing that the so-called "Hubbard-I approach" [^0) leads to 
a selfenergy being identical to that of the W = limit. It thus coincides with 
that of the SDA for B k -„ = 0. Comparing the results of SDA and Hubbard-I 
therefore helps to understand the implications due to B k -a- 

Two severe shortcomings of the SDA are also very obvious. The first is the 
neglect of quasiparticle damping by the delta- function ansatz (34) which leads 
to a real selfenergy £§® A (E). The second concerns the weak-coupling regime, 
which is surely violated by the strong-coupling theory SDA. It is not to expect 
that Fermi liquid behaviour is correctly described by the SDA. These two points 
shall be attacked and eliminated by the subsequent approaches. But let us first 
inspect what the SDA tells about bandferromagnetism in the Hubbard model. 

For a full SDA-solution we have to determine the "band correction" Bk_ a (9). 
The local part B- a does not need a special treatment because of its direct con- 
nection (11) to the spectral density itself. The fc-dependent part Fk_ a appears 
a bit more complicated. Assuming translational invariance and next neighbour 
hopping, only, the fe-dependence can be separated ]l7|,^8[: 



(n_ ff )(l - (n_ <7 ))F fc _ . = (e(fe) - 


-To)J2 Fi -l 


(36) 


FrJ = (m- a rij- a ) - {n-a) 2 




(37) 


r — a = \ c j- a c j a c i-o- c io-) 




(38) 


* ' ~<J = \ C jo C J-° C i-a C i<?) 




(39) 



Ferromagnetism in the Hubbard Model 9 

i,j are numbering next neighbours. The method used to determine the higher 
correlations F_„ shall be exemplified for F_l- First we rewrite F_„ as 

F -i = E S fl ^J-aC+ +A _ a C ] + Aa ) (40) 



where the index A corresponds to the lattice vector which connects two neigh- 

(3) 

bouring sites Ri and Rj . Because of translational invariance F_J. does not de- 
pend on the explicit value of A. We now introduce a "higher" spectral density 

(3) 

Sj s J(E) as the (fc, ^-dependent Fourier transform of 

sfAt,t>) = ^([(c^c+ +A _ a c 3+Aa )(t),ci(t')} + ) (41) 

The spectral theorem yields: 

oo 

F -* = M ^ / dEf - {E)S ^ {E /i} (42) 

k -oo 

According to the definition (41) the poles of Sj^J(E) must belong to the single- 
electron excitations of the Hubbard system. From the spectral representation 
of SfrJ (E) and by comparison with the respective representation of the single- 
electron spectral density Sk a (E) again a two pole ansatz appears to be consis- 
tent, 

2 

S ktr (E) = n^2& jtr (,k)6(E + n-E jtr (k)) (43) 

3=1 

where the quasiparticle energies Ej a (k) are the same as those in (34) so that the 
spectral weights a.j a are the only unknown parameters. They are fixed by the 
first two spectral moments of Sj^J (E) leading via (42) then to an (approximate) 
expression for F_J.. 

In an analogous manner the correlation terms F_^ and -F_ CT are determined 
by two-pole ansatzes of properly chosen "higher" spectral densities S^ (17| . 
We note in passing that the same procedure can of course also be applied to the 
term (ni a c^_ a Cj- a ) in the "band shift" B- a (10) yielding then the exact result 
(11). 

Fig. 2 demonstrates for an sc lattice that the just-described SDA allows for 
bandferromagnetism in the Hubbard model. The zeros of the inverse param- 
agnetic susceptibility \~ l as a function of the band occupation n indicate in- 
stabilities of the paramagnetic state towards ferromagnetism (T = OK). The 
instabilities appear as soon as U exceeds the critical value Uq ~ 4W. It is a 
special feature of the SDA JlJ],[l9| , maybe even of the Hubbard model itself, that 
there appear for U > Uc two ferromagnetic solutions, i.e. two zeros of X ■ The 
first solution sets in at nc > 0.34 where the actual value only slightly depends on 
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0.4 0.6 
n 

Fig. 2. Inverse paramagnetic static susceptibility x _1 f° r the sc lattice as a function 
of the band occupation n for various values of U. (a) System with the full fc-dependent 
selfenergy. (b) System with a local selfenergy (Fk- a = 0). The inset shows the magnetic 
moment as a function of n for U/W = 5. 



U running into saturation (m=n) for n > 0.68 (see inset Fig. 2). The second so- 
lution appears for higher band occupation, but does never reach saturation and 
is always less stable than the other solution. Two important aspects of the SDA 
solution in Fig. 2 should be stressed. The hrst concerns the comparison with the 
corresponding results of the Hubbard-I approach |I(J , which does not allow ferro- 
magnetism on the sc lattice. On the other hand, the Hubbard-I selfenergy differs 
from (35) "only" by neglecting the "band correction" B k - a - Consequently, the 
stronger magnetic stability in the SDA must be due to B k - a . The second remark 
aims at the non-locality (fc-dependence) of the electronic selfenergy provoked by 
the "bandwidth correction" F k - a (36) in B k -„. Part (b) of Fig. 2 demonstrates 
its importance. Switching off this term leads to a dramatic further increase of 
the critical coupling ^ from 4 to 14. More detailed studies |y|, however, show 
that the influence of F^-a on magnetic stability strongly depends on the lattice 
type. Increasing coordination number (sc — > bcc — ► fee) drastically diminishes 
the importance of the non-locality leaving the local spin-dependent band shift 
B- a as the pushing mechanism for ferromagnetic stability. 

It can be shown |l7| that the SDA gives a qualitatively convincing picture of 
bandferromagnetism . The main message is that the spin-dependent "band-shift" 
B- a and the lattice structure are the most important ingredients. However, we 
should not forget the disadvantages of the method. One of them is the neglect of 
quasiparticle damping, which is assumed to destabilize the collective spin order. 



3.2 Quasiparticle Damping 

More or less by construction the SDA selfenergy (35) is real except for a single 
5-peak in Im£ ka (E + i0+) at the pole E = T + B k ^„ + U(l - (n- a )) - /i. It 
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has, however, no influence since it falls into the Hubbard gap. The quasiparticles 
are stable. We are therefore looking for an approach that retains the obvious 
advantages of the SDA but improves it by a reasonable inclusion of quasiparticlc 
lifetime effects. 

Starting point is an alloy-analogy for the Hubbard model which traces back 
to Hubbard himself p3] . The idea is to consider the propagation of a a electron 
through the lattice with the —a electrons being "frozen" at their lattice sites 
and randomly distributed over the crystal. When the a electron enters a lattice 
site it can encounter two different situations. It can meet a — a electron or not. 
That can be interpreted as a hopping through a fictitious binary alloy, the two 
constituents of which are characterized by atomic levels E\ ai E^a and by concen- 
trations xia-,X2a- The "coherent potential approximation" (CPA) is a standard 
method to perform the configurational average over the alloy |2l| . The resulting 
a selfenergy obeys the following implicit equation |ll| , 

E pa — £(j (E) — T 



p=\ 



l-±G a {E){E pa -E„(E)-T Q ) 



G a (E) — jj^2Gka(E) is the local propagator. CPA is a single-site approxi- 

k 

mation, the resulting selfenergy therefore wave vector independent. To proceed 
one has to specify the two-component alloy. The natural way ("conventional al- 
loy analogy" CAA) would be to use the zero-bandwidth results of the Hubbard 
model: E\ a = Tq,E2o- = Tq + U,xi,j = 1 — (n- a ) = 1 — X2a- The resulting 
■^a (E) includes quasiparticle damping but excludes any spontaneous ferro- 
magnetism |22|,^3[ . Does quasiparticle damping really kill any spontaneous mo- 
ment order? This is a serious question because it has been shown |24j that for 
infinite lattice dimensions CPA is an exact (!) treatment of the alloy problem. 
On the other hand, the CAA selfenergy disagrees with the strong as well as 
weak coupling behaviour (24), (28). This discrepancy can only be removed by 
the conclusion that the zero-bandwidth limit constitutes the wrong alloy analogy 

The fictitious binary alloy is indeed by no means predetermined. We consider 
E\ t 2a and Xi^a at first as free parameters and fix them via the high-energy 
expansions (21) and (24), which we insert into the CPA equation (44). Then we 
expand eq. (44) with respect to powers of 4 . Comparing coefficients we can use 
the following set of equations to determine an "optimum alloy analogy" : 



/ %pa J- 




P =i 




2 




/ •Ep(7\l- J p<j 


-T Q ) = U(n_ a ) 


p=l 




2 




/ Xpa \t^pa 


-To) 2 = U 2 {n- a ) 



p=l 
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2 

5> pff (£ pCT - T ) 3 = U 3 (n_ a ) + U 2 B_ a (n_ a ){l - (»_„)) (45) 

P =i 

Deriving from these equations E pa and x po - automatically guarantees the correct 
strong-coupling behaviour. One finds (p=l,2) ("modified alloy analogy" MAA). 



To + 2 P + B -v + (-l)V(^ + B-a) 2 - 4(n_ CT )S_ ;7 ] 

= (E^ DA (k)) e(k)=To = / P (T , J7, (n_ CT ), B_ CT ) (46) 



maa _ -B-g + T + £7(1 - (n_ g )) - E^ AA _ _ AMj4 , . 

^ AA = (c4? A ) e(fe)=To = , 9p (T , C7, («_,), S_ CT ) 

Surprisingly, the energies and weights coincide exactly with the corresponding 
SDA entities if e(k) is simply replaced by T . Because of the single-site aspect of 
the CPA pl|jLll the "band-correction" B k -a is restricted here to its local part 
-B-o-, the decisive band shift. Inserting (46) and (47) into (44) yields the MAA 
selfenergy for all E which exhibits some remarkable features pa : 

(1) As a CPA result the MAA includes quasiparticle damping (ImS^ IAA (E) =/= 
0), and that without giving up the advantages of the SDA. 

(2) The expectation values (n— a ) and -B_ CT are to be determined selfconsistently. 
In principle, they take care for a carrier concentration-, temperature- and 
spin-dependence of the atomic data (46), (47) of the alloy constituents. Fur- 
thermore, and that is quite an important aspect, B_ a brings into play in a 
certain sense the itineracy of the — a electrons ("correlated" electron hop- 
ping), completely neglected in the CAA. 

(3) Strong-coupling and high-energy behaviour are correctly reproduced, more 
or less by construction. 

(4) The general CPA theory |2l]] comes for the so-called split-band regime (here, 
U ^S> W) to the conclusion that the spectral density Sk a (E) should consist 
of two separated peaks with centres at 

TCf A = E pa + x pa (e(k)-T ); p = 1,2 (48) 

Inserting (46) and (47) for U 3> W yields exactly the Harris-Lange results 
(6-8), a further strong confirmation of the MAA. 

(5) Contrary to the CAA the MAA allows for spontaneous bandferromagnetism. 

For a typical example the possibilities of the MAA are demonstrated in Fig. 3, 
which shows Sk (7 {E) for strongly correlated electrons {Uf = 5) on an fee lattice. 
For less than half-filled bands (n < 1) the system is paramagnetic. The band oc- 
cupation n=1.6 used in Fig. 3 allows for a spontaneous collective order provided 
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Fig. 3. Spectral density as a function of energy for an fee lattice calculated within the 
MAA, (a) T = WOK, (b) T = 615K for different fe-vectors equidistant along the (001)- 
direction of the 1. Brillouin zone. Further parameters: n — 1.6, U — 20eV, W = 
4eV. The vertical line indicates the position of the chemical potential. 



U exceeds a critcal value. Two types of splitting occur. At first the spectral den- 
sity consists for each A;- vector of a high-energy and a low-energy peak separated 
by an energy of order U. The finite widths of the peaks are due to quasiparticle 
damping and obviously energy-, wave- vector-, spin- and temperature-dependent. 
The spectral weight (area) of the low-energy peak scales with the probability that 
the propagating (ft, cr)-electron in the more than half- filled band enters an empty 
lattice site, the weight of the upper peak that it meets anywhere a — a electron. 
This splitting is general and not at all bound to ferromagnetism. It demonstrates 
the correct strong-coupling behaviour (6,7). Ferromagnetism appears when the 
two spectral density peaks show a spin asymmetry. For the low temperature 
example T = 100K (Fig. 3a) the system is close to saturation (m = 2 — n), 
i.e. the up-spin states are almost fully occupied. A down-spin electron cannot 
avoid to meet an up-spin electron at every lattice site being therefore forced to 
perform a Coulomb interaction. Consequently the low-energy peak of Ski(E) 
disappears. At higher temperatures (Fig. 3b) the peak reappears because of a 
partial demagnetisation, i.e. a finite density of holes in the up-spin spectrum. 
At low temperatures the high-energy peaks of Sk[(E) are very sharp indicat- 
ing long-living quasiparticles. A J.-hole has no chance to meet an j-hole to be 
scattered. 

An interesting fe-dependence of the exchange splitting is observed. At the 
branch-top (X-point) a "normal" splitting apppears, i.e. the down-spin peak is 
above the up-spin peak. At the bottom (/ n -point), however, the t-peak is higher 
in energy than the |-pcak ("inverse exchange splitting"). The quasiparticle dis- 
persions of the two spin-parts are crossing as a function of k. This behaviour 
is a result of two competing correlation effects. The one is the spin-dependent 
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exchange shift of the centres of gravity of the quasiparticle subbands, the other 
a spin-dependent band narrowing. The latter may overcompensate the first. In 
any case the exchange splitting exhibits a strong wave-vector dependence, even 
with a sign change. 

We can conclude that the MAA improves the SDA in a systematic manner 
by including quasiparticle damping. As will be discussed in Sect. 4, the main 
modification consists in a substantial destabilisation of ferromagnetism in the 
Hubbard model. One reason is the finite overlap of the spectral density peaks, 
in the SDA because of (34) excluded, which weakens the ferromagnetic solution 
of the selfconsistently evaluated model theory. 

The second shortcoming of the SDA, the incorrect weak-coupling behaviour, 
persists in the MAA. In the following we have to inspect its consequences for 
the strong-coupling phenomenon ferromagnetism. 



3.3 Weak-coupling behaviour 

Neither the SDA selfenergy (35) nor the MAA selfenergy does fulfill the weak- 
coupling expansion (28). We are now looking for an approach which interpolates 
reasonably between the strong- and weak-coupling regimes. It should extend the 
preceding approaches SDA and MAA by inclusion of weak coupling and Fermi 
liquid behaviour without giving up the convincing essentials from the SDA and 
MAA: spin-dependent band shift B- a and quasiparticle damping. Following an 
idea of Kajueter and Kotliar [[L5| we start with a selfenergy ansatz |2{|, which 
we call the "modified pertubation theory" (MPT): 

Ska{E) = U{n- a ) + 2 — t=- (49) 

The second order contribution A ka (E) is defined in eq. (29). We use in A ka (E) 
for the spectral densities the Hartree-Fock version (30) with the selfconsistency 
condition (31). By construction the selfenergy (49) is correct in the weak-coupling 
regime up to order U 2 . In order to fulfill simultaneously the strong-coupling 
regime we expand (49) in powers of the inverse energy 1/E and compare this 
with the exact selfenergy expansion (24). The first two terms are automatically 
fulfilled, the third and the fourth fit the coefficients a klJ , b ka @1- For the results 
presented in the next section we have implemented the MPT (49) in a "dynamical 
mean field theory" (DMFT) procedure [ J26|]27| . The main assumption is then 
a local selfenergy S ka (E) — > E a (E) , exact in infinite lattice dimensions Q. 
Furthermore, we exploit the fact that the Hubbard problem can be mapped under 
such conditions on the single-impurity Anderson model (SIAM) (see contribution 
of D. Vollhardt ) as long as a special selfconsistency condition is fulfilled pq| . 
The idea of the MPT is then applied to the simpler SIAM problem. It can be 
shown, however, that the direct application of the MPT-concept to the Hubbard 
model gives almost the same results. 

The MPT fulfills a maximum number of limiting cases: The weak-coupling 
behaviour is correct up to C/ 2 -terms, for all bandoccupations n Fermi liquid 
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Fig. 4. Curie temperature as a function of band "filling" (hole concentration) for an 
fcc-type d = oo-lattice (51). The points with error bars are the Quantum-Monte-Carlo 
results from ref. [3]. SDA, MAA and MPT are explained in the text. 



properties are recovered. The Luttinger theorem, 



SuA0) = S^(0) 



l-j. = m|u=o + E<r(0) 



(50) 



is fulfilled for all bandoccupations n, at least if not too far from half-filling 
(n = 1). For low temperatures a Kondo resonance appears at the chemical po- 
tential /i. The zero-bandwidth limit (W — ► 0) comes out exact for all n as well 
as the strong-coupling behaviour gathered in eqs. (6)-(8). Eventually bandferro- 
magnetism is possible within the framework of MPT. According to the number 
of correctly reproduced limiting cases the MPT appears to be an optimum ap- 
proach to the many-body problem of the Hubbard model. 

4 Discussion 



Let us compare the analytical approaches presented in the preceding Section with 
respect to the most important magnetic key-quantity, the Curie temperature Tq. 
For comparison with the numerically essentially exact Quantum Monte Carlo 
results of Ulmke fl we have evaluated the three theories for an fcc-d = oo type 
lattice described by the Bloch density of states: 



Po(E) 



exp(-!(l-^)) 



PWttCL-^) 



(51) 



The energy unit is chosen to t* = t\/2d{d — 1) = 1. The density of states (51) 
is strongly asymetric with a square root divergency at the upper edge. 

All the presented methods predict that ferromagnetism exists only for more 
than half-filled energy bands. In Fig. 4 bandfilling therefore means the hole- 
density and the DOS is that for holes following from (51) by t* — » — t*. 
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Fig. 5. Left: Quasiparticle density of states as function of energy for the fee-type lattice 
(51). Solid lines for a —"\, dotted lines for a =J,. QMC results from [3]. Temperature 
T is chosen so that for all presented methods a magnetization m = 0.4 is found. 
Right: For the same parameters as in the left part the magnetisation m and the inverse 
paramagnetic susceptibility x _1 as function of the reduced temperature T/Tc- Bottom 
figure: QMC results, dashed line for T/T c < 1 : 5 = 1/2-Brillouin function. 



The role of the spin-dependent band shift _B_ CT (10) for the magnetic stabil- 
ity becomes evident by comparison of the three methods with their B- a = 
counterparts. That is for SDA the Hubbard I-solution |l(| which yields ferro- 
magnetism only for very asymmetric DOS. In the case of (51) there appears 
ferromagnetism in a strongly restricted region of very low hole densities. The 
band shift B^ a in the SDA obviously leads to a drastic enhancement of ferro- 
magnetic stability. The counterpart of MAA is the conventional alloy analogy 
(CAA) that does not allow for any n ferromagnetic order |g^,g3|. Finally, the 
MPT-counterpart, that neglects the bandshift, is the Kajueter-Kotliar approach 
[ [15| . It indeed exhibits ferromagnetism but with substantially lower Curie tem- 
peratures To 

The influence of quasiparticle damping can best be judged by comparing the 
results of SDA and MAA. The inclusion of lifetime effects in MAA presses the 
Curie temperature to half the SDA-values. The ferromagnetic coupling strength 
is therefore substantially weakened by quasiparticle damping. The incorporation 
of the correct weak-coupling behaviour by MPT leads to a further Tc-reduction 
possibly because of the screening tendency with respect to the effective mag- 
netic moments manifesting itself in the appearance of a Kondo resonance. The 
MPT-results for Tc(n) are closest to the essentially exact numerical results of 
Ulmke H , which have been read off from the zeros of the inverse paramagnetic 
susceptibility. 
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The similarities and the differences of the various methods come out by the 
quasiparticlc density of states in Fig. 5a. All theories show the splitting into 
the two Hubbard bands and the additional exchange splitting which takes care 
for the spin asymmetry. As a consequence of the neglect of damping effects the 
DOS-structures are rather sharp in the SDA. They appear smoother in MAA. 
The correct low-energy behaviour of the MPT leads to the best approach to the 
Q-DOS found by "Quantum Monte Carlo" (QMC)-evaluation of a "Dynamical 
Mean Field Theory" . Note that for a reasonable comparison of the Q-DOS we 
have chosen in Fig. 5a for the various theories temperatures which lead to the 
same magnetization m = 0.4. 

Fig. 5b demonstrates impressively the qualitative equivalence of the presented 
methods when we plot the spontaneous magnetization and the paramagnetic 
susceptibility as function of the reduced temperature T/Tc- There are indeed 
strong similarities. All magnetization curves are Brillouin function- like with full 
polarization at T = 0. The only exception is MPT, which shows up a slight 
deviation from saturation at T = 0. For the chosen parameter set all methods 
predict a second order transition at Tq ■ The paramagnetic susceptibility follows 
in all theories a ideal Curie- Weifi behaviour: \ = C(T— 0)~ l . The paramagnetic 
Curie temperature equals in any case the Curie temperature Tq and even the 
Curie constants C are all very similar: C =0.42 (SDA), 0.52 (MAA), 0.57 (MPT), 
0.47 (QMC). One important consequence of this striking qualitative equivalence 
and correctness of SDA, MAA and MPT compared to QMC is that even the 
rather simple SDA can be used to describe the magnetism of more complicated 
structures (films, surfaces, multilayer, real substances). 

We conclude that ferromagnetism does exist in the Hubbard model depending 
on lattice structure, band occupation, Coulomb coupling and temperature. By a 
critical comparison of three different analytical approaches we could demonstrate 
the essentials for ferromagnetic stability, in a positive sense the spin-dependent 
band shift B- a , which represents the correlated hopping of opposite spin elec- 
trons, and in a negative sense the finite lifetime of quasiparticles, which results in 
a drastic lowering of the ferromagnetic coupling strength. The correct inclusion 
of Fermi liquid properties eventually leads to the most convincing description of 
ferromagnetism in the Hubbard model. 
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